www.gusucode.com > 小波分析理论与Matlab 7实现 > 小波分析理论与Matlab 7实现/code/Ch16/examp16_23.m
%装载并显示原始图像 load bust; %加入含噪 init=2055615866; randn('seed',init); X1=X+20*randn(size(X)); subplot(221);image(X1); %colormap(map); title('原始图像');axis square; %用小波db4对图像X进行一层小波包分解 T=wpdec2(X1,1,'db4'); %重构图像近似部分 A=wprcoef(T,[1 0]); subplot(222);image(A); title('图像的近似部分'); axis square; %边缘检测 %%原图像的边缘检测 BW1 = edge(A,'sobel'); subplot(223);imshow(BW1); title('原图像的边缘'); axis square; %图像近似部分的边缘检测 BW2= edge(X1,'sobel'); subplot(224);imshow(BW2); title('图像近似部分的边缘'); axis square;